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We analyze the CP asymmetries of B — > and 1? — > rf''K modes in the QCD improved fac- 
torization framework. For our calculation we use the phenomenological parameters predetermined 
from the global fit for the available B — ► PP and VP modes (without the quark-level subpro- 
cess b — > sss). We show that the large negative sin(20i)0_K- a and the large branching ratio for 



> 

qq ! B^ — ► n' can be simultaneously explained in the context of supersymmetry (SUSY). The R- 

m 

psj ' parity conserving SUGRA models are used and their parameter space is constrained with the 

7-H ! 

C*~) ■ observed dark matter relic density along with other experimental constraints. The R-parity vio- 

o. 

lating SUSY models are also used to show that they can provide solutions. We calculate the CP 



asymmetries for different B^C ) — > <f)K^^ and B^ ) — * rf^K^-^ modes and show that the SUSY 
model predictions are consistent with the available experimental data. 
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Recent results from Belle and BaBar on B decays provide an opportunity to test the 
source of CP violation. The standard model (SM) source for CP violation arises from the 
CKM matrix which has only one phase. The SM predictions for CP asymmetry for different 
B decay modes are being tested at the B factories and discrepancy seems to be emerging. 
For example, the modes B — > <pK s and B — > J /^>K S are uniquely clean in their theoretical 
interpretations. Among these decay modes, B — > (f)K s occurs only at one loop level in the 
SM and hence is a very promising mode to see the effects of new physics. In the SM, it is 
predicted that the CP asymmetries of B — > <j)K s and B — > J/^K S should measure the same 
sin(20i) with negligible 0(\ 2 ) difference The Belle and BaBar experiments measure: Q] 

sin(20i) J/ * iC - J = 0.734 ± 0.055, (1) 

and 0,Q 

sin(20!) 0Ks = -0.96 ± 0.51°;°? (Belle), 0.45 ± 0.43 ± 0.07 (BaBar). (2) 

The world average shows a 2.7cr disagreement between sin(20i)^ s and sin(20i) j/^k 3 - The 
sin(20i) j/^k s being a tree level process is in excellent agreement with the SM theoretical 
prediction, sin(20 1 ) SM = 0.715±g;gti B 

Experimental data is also available for the decay modes B^ ' — > ^(')j^±(o) w hich involve 
the quark- level subprocess b — > sss as in B — > <pK s . The world average value of the measured 
branching ratio (BR) is B(B ± -> rj'K^ = (77.6 ± 4.6) x 10~ 6 

BBS 

, which is larger than 

the predicted SM value. The results also exist for sin(20i)^x s = 0.33 ± 0.34 9l lld| and the 
CP rate asymmetries Acp for different B^ — > (j)K ± and B ± — > 7]^K ± modes 

In this letter we try to find a consistent explanation for all the observed data in charmless 
hadronic B — ► PP and B — > VP decays [-P( V) denotes a pseudoscalar (vector) meson] in the 
framework of QCD factorization (QCDF). We calculate the BRs and the CP asymmetries 
for the decay processes B±^ _> ^±(0) and B±^ -> rj^K^ in the SM and its SUSY 
extensions with R-parity conservation (SUGRA models), and with R-parity violation. The 
required QCDF input parameters for the calculation are determined by using a global fit of 
all possible B decay modes without the subprocess b —>■ sss, since this process may involve 
new physics. 

Previously, new interactions were invoked to exp lain the large BRs of B ±( - ^ — > rj'K^^ 



and the large negative sin(20!)0x s [I2ul3tll4ll5lll6ul7l ]. However, attempts of simultaneous 



explanation for the observed data of these decay modes were made only by using the naive 



actorization technique 



13, 



151] - We adopt the newlv developed OCD improved factorization 
3| for the calculation in this work. This approach allows us to include the possible non- 
factorizable contributions. In the heavy quark limit » Aqcdj the hadronic matrix 
element for B — > M\M 2 due to a particular operator Oj can be written in the QCDF as 
follows: 



(M 1 M 2 \O i \B) = {M l M 2 \O i \B)^ F 



Aqcd 



(3) 



where NF denotes the naive factorization. The second and third term in the square bracket 
represent the radiative corrections in a s and the power corrections in Aqcd/^6- The decay 
amplitudes for B — > M\M 2 can be expressed as 

A(B -> M X M 2 ) = A f (B — > M X M 2 ) + A a (B -> M X M 2 ) , (4) 



where 



G 10 

A f (B — > MiM 2 ) = -| £ E^*^M a i( M l M 2l°il B >NF 

;,c i=i 



Z^ • P" ' pg 

V * p=u,i 
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^(5 - M X M 2 ) = ^ / b / Mi /m 2 E E W ^ • (5) 

V Z p=u,c i=l 

Here A^(B — ► M\M 2 ) includes vertex corrections, penguin corrections, and hard spectator 
scattering contributions which are absorbed into the QCD coefficients a^, and A a {B — > 
Mi Mo) includes weak annihilation contributions which are absorbed into the parameter 

end-point divergent integrals are parametrized as [18|] 



/ox ^ ' A h 

where Xh and denote the hard spectator scattering contribution and the annihilation 
contribution, respectively. Here the phases 4>h,a are arbitrary, 0° < 4>h,a < 360°, and the 
parameter Ph,a < 1> an d the scale A h = 0.5 GeV assumed phenomenologically [18]. In 
principle, the parameters Ph,a an d <Ph,a for B —>■ PP decays can be different from those for 
B — > VP decays. 

We use the global analysis to determine the QCD parameters Xh and Xa Q|- I n the 
global analysis we exclude the decay modes whose (dominant) internal quark-level process 



is b — > sss: for example, B — > 0if and B — > rj^'M, where M denotes a light meson, such as 
7T, -R', p, if*. The reason for such a global analysis is due to the fact that the b —>■ sss mode 
may require the existence of new physics. We use twelve B — > PP and VP decay modes, 
including B — > nn, irK, pn, pK, un, uK, to determine the QCD parameters, pYh^ and 
4>a P h PV ■ The parameters are distinguished by their superscripts P and V which denote final 
state mesons. Now, if pa,h and 4>a,h are large, the effects of Xa,h can be large and give rise 
to a large annihilation (and/or hard spectator scattering) contribution, which indicates the 
theory becomes less reliable and suspicious due to that large non-perturbative contribution. 
We find that it is possible to obtain a good global fit (Xmm = 7-5) with large Xa,h effects. In 
this scenario, the BR for B + — > rj'K + is 74.7 x 10~ 6 , which saturates the experimental limit 



with the very large Xa,h effects jl7| and the BR for B + — > <pK + is small: 4.0 x 10~ 6 . It is 
possible to obtain successful fits to all B — > <pK and B —>■ rj'K data in our SUGRA models. 
But, since the effects of Xa,h are large, we will just comment on the fits in this scenario in 
our result sections. 

In Ref. jlfl ]. we generated a fit (with xLin = 18-3) with the (relatively) small Xa and 
X H effects. The corresponding theoretical inputs for this fit are as follows (we will use these 
values in this work): 

A = 0.2198, A = 0.868, 3 = 86.8°, |V u6 | = 3.35 x lO" 3 , 

p = 2.1 GeV, m s (2 GeV) = 85 MeV, f B = 220 MeV, 
F Bn = Q 2A9 ^ R ^ = 1; a b p = 31; 

p7 = 0, P r = 0.5, p p H p = l, p v / = 0-746, 

<PY = -6°, ^ = ^ = 180°. (7) 

Note that in this case the effect of the weak annihilation parameter Xa is relatively small 
(i.e., pY = and p\ p = 0.5), and the effect of the hard spectator scattering parameter X H 
is also very small, because p|J p = 1, p^ p = 0.746, and 0ff P = 0^- p = 180° so that the terms 
1 and pne^ H in Xjj cancel each other in Eq. ©. Based on the above inputs, the BRs and 
CP asymmetries for B — ► <pK and B — > rj'K are predicted and shown in Table I. We see that 
the predicted B(B + ^ — > r)'K + ^) are smaller than the measured values and the predicted 
sin(20i)^x s is also very different from the world average. The predicted B(B + — > i]K + ) is 
smaller than the experimental data as well. 

In the following two sections, we will discuss the CP asymmetries and BRs of B — > (j)K 



TABLE I: The branching ratios (in unit of 10 ) and CP asymmetries of B - 
decays are shown. The values in the parenthesis are experimental numbers P, U, 0, 7 



K and B -^_n^K 
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I20I ]. Here the inputs for (He lit. wit 1 1 the small eliecls of A* 1 and A// are used. 'I lie lirsl row of the 
last column shows the sin(2</>i) for B —* </>if s and the last row in the same column shows the value 
rfK 8 . 
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^CP 


Decay mode 


BR 


^CP 


sin (20i ) 


B + -> 


7.3 





B° -> <j)K° 


6.7 


0.01 


(j)K s : 0.68 




(9.2 ±1.0) 


(0.03 ± 0.07) 




(7.7 ± 1.1) 


(0.19 ±0.68) 


(Eq. H 


S+ -> 7?'Er + 


51 


0.01 


B° r/if 


46.8 


0.016 


r/'K s : 0.57 




(77.6 ±4.6) 


(0.02 ± 0.04) 




(60.6 ± 7) 


(0.8 ±0.18) 


(0.33 ± 0.34) 




1.9 


-0.16 




1.7 


-0.16 






(3.1 ±0.7) 


(-0.32 ±0.20) 
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(-) 





and B — > rj'K modes in the context of SUSY models. 
[1] R-parity violating SUSY case 

The R-parity violating (RPV) part of the superpotential of the minimal supersymmetric 
standard model can have the following terms 



W RPV = KiLiH 2 + UjkLiLjEZ + l' ijk LiQjDl + I'^Dp 



c 

k i 



(8) 



where E iy Ui and Di are respectively the z-th type of lepton, up-quark and down-quark 
singlet superfields, Li and Qi are the SU(2)^ doublet lepton and quark superfields, and H 2 
is the Higgs doublet with the appropriate hypercharge. 

For our purpose, we will assume only /'—type couplings to be present. Then, the effective 
Hamiltonian for charmless hadronic B decay can be written as [21], 



H eff( 



^jkn 



u jkn 



-I- rl L 

' u jkn 



Here the coefficients dj^ and u%„ are defined as 



(9) 



jkn 



l jkn 

3 // 7/* 

E tj ijk i in'i 
1 8ml - ' 

t=l v iL 



"'jkn 



1=1 



5m 



2 ' 



u jkn 



3 ]i ]>* 



n 



2) 



(10) 



i=l ^""e iL 

where a and /3 are color indices and 7^ L = 7 M (1 ± 75). The leading order QCD correction 



to t 



21 



lis operator is given by a scaling factor / ~ 2 for 777.5 — 200 GeV. We refer to Refs. 



221 ] for the relevant notations. 



The RPV SUSY part of the decay amplitude of B — > <pK is given by 



A 



RPV 

<j>K 



(rt L _|_ fi 

I u 222 1 u : 



222 



where the coefficient a is expressed as 

C F a s 



1 

Ct = -r-r- 

It has been noticed 
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(ii) 



;i2) 



22j that the RPV part of the decay amplitude for B — > r/'iT, ^4S^, 



is proportional to (g? 222 — <if 22 ), while the RPV part of the decay amplitude for B — > <pK, 
, is proportional to (d 222 + <if 22 ). It has been also pointed out [22j that the opposite 
relative sign between <i 222 and in the modes B — > r/'iT and B — > appears due to 
the different parity in the final state mesons 7/ and 0, and this different combination of 
(rf 222 — d 222 ) and ( ^222 + ^222) m these modes plays an important role to explain both the 
large BRs for B — > rj'K and the large negative value of sin(20i)^ s at the same time. 
We define the new coupling terms (i 222 and d 222 as follows: 



^222 ^ |A„"qoA 



«32' v «22l 



J0L 



G? 222 OC |/X^9q/\ 



i22' v i23l 



(13) 



where 9i and 9r denote new weak phases of the product of new couplings A- 32 A^ 22 and 
^22^*23; respectively, as defined by A 332 A 3 * 22 = |A 332 A 3 * 22 |e^ and A 322 A 3 * 23 = |A 322 A 3 * 23 |e iefl . 

We consider two different follows. 

Case (a): d 222 7^ and d 222 = 0, i.e., 



IA322I - 0.077 , |A 332 



0.077 , |A 323 | 
6 L = 1.5 , msusv = 200 GeV. 



(14) 
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Our results are summarized in Table II. We use the negative d 222 and the scaling factor 
We find that sin(20i)0x s can be brought down to at most. The BRs, B(B + — > r)K + ) and 
B(B° — > T]K°), are larger compared to the experimental values, but the experiments (Belle, 
BaBar and CLEO) are not quite in agreement and the experimental errors are also large 



TABLE II: Case (a): the branching ratios (in unit of 10 6 ) and CP asymmetries of B — ► <pK and 
B — > rj^K decays are calculated in the framework of R-parity violating SUSY. 
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*4cp 


Decay mode 


BR 


A C P 


sin (2^i ) 


B+ -> 




8.9 


-0.17 


B° -» (j)K° 


8.2 


-0.18 


<pK s : -0.03 






72.0 


0.16 


B° -> r/if 


66.0 


0.16 


rfK a : -0.2 


5+ -» 


7]K+ 


10.5 


0.25 




9.7 


0.25 





TABLE III: Case (b): the branching ratios (in unit of 10 6 ) and CP asymmetries of B — > <f>K 
and -B — > 77W X decays are calculated in the framework of R-parity violating SUSY. 
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BR 


Ac? 


Decay 


mode 


BR 


A C P 


sin (20i ) 


B + -> 


10.2 


-0.05 


B ^ 


<pK Q 


9.5 


-0.04 


(pK s : -0.61 


5+ -> 7?'iC + 


74.0 


0.11 


B ^ 


V 'K° 


67.7 


0.11 


7]'K S : 0.48 


5+ 7^+ 


6.7 


0.06 


B ^ 


7]K° 


6.1 


0.06 





for these modes. In Table II, we used 5' = 0, but if we use 5' = 30°, the sm(2(j)i)^ Ks can be 
larger negative: —0.2, after satisfying all the constraints [especially the B(B + — > r]K + ) and 
the sin(20i)^J. 

Case (b): d% 22 ^ dg 2 ^ 0, i.e., 

|A' 322 | = 0.076 , |A 332 | = 0.076 , |A' 323 | = 0.064 , 

9 L = 1.32 , 6 R = -1.29 , msusv = 200 GeV. (15) 

Our results are summarized in Table III. We find that sin(2<f)i) ( f > x s can be large negative. 
In addition to the parameters given in Eq. (fTHjl . we also used the additional strong phase 
5' = 30°, which can arise from the power contributions of h.QCD/ m b neglected in the QCDF 
scheme, and whose size can be in principle comparable to the strong phase arising from the 
radiative corrections of 0(a s ). If 5' = is used, we obtain sin(20i)^ s = —0.2. 

In the case of the large X^h effects with xLin = 7-5 where the BR of B + — > r]'K + is large, 
we can use the R-parity violating SUSY couplings to raise the BR of B + —>■ <j)K + (which is 
small, 4.0 x 10~ 6 to begin with). It is possible to raise B(B + -> <j)K + ) to (8 - 9) x 10" 6 , but 
its CP rate asymmetry Acp{B + — > <pK + ) is large ~ —0.4 and sin (20i)^ Xs can be brought 
down to at most —0.16. 
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The RPV terms can arise in the context of SO(10) models which explain the small 
neutrino mass and has an intermediate breaking scale where B — L symmetry gets broken 
by (16 + 16) Higgs. These additional Higgs form operators like 16#16 m 16 m 16 m /M pi (16 m 
contains matter fields) and generate the RPV terms j^ . 

[2] R-parity conserving SUSY case 

As an example of the R-parity conserving (RPC) SUSY case, we will consider the su- 
pergravity (SUGRA) model with the simplest possible non-universal soft terms which is the 
simplest extension of the minimal SUGRA (mSUGRA) model. In this model the lightest 
SUSY particle is stable and this particle can explain the dark matter content of the universe. 
The recent WMAP result provides j^J: 

{l CDM h 2 = 0.11261°;°°*, (16) 

and we implement 2a bound in our calculation. 

In the SUGRA model, the superpotential and soft SUSY breaking terms at the grand 
unified theory (GUT) scale are given by 



,2| i |2 



W = Y U QH 2 U + Y D QH 1 D + Y L LH 1 E + nH x H 2 , 

'1 



V m a X a X a + Bfj,HiH 2 



£soft = ~X! m i 

i 

+ (A U QH 2 U + A D QH 1 D + A L LH 1 E) + H.c. , (17) 



where E, U and D are respectively the lepton, up-quark and down-quark singlet superfields, 
L and Q are the SU(2)x, doublet lepton and quark superfields, and H X2 are the Higgs dou- 
blets, (pi and X a denote all the scalar fields and gaugino fields, respectively. The parameters 
in the mSUGRA model, a universal scalar mass mo, a universal gaugino mass mm, and the 
universal trilinear coupling A terms are introduced at the GUT scale: 

m\ = ml m a = m 1/2 , A U ' D ' L — A Y U > D,L , (18) 

where Y U,D,L are the diagonalized 3x3 Yukawa matrices. In this model, there are four 
free parameters, m , mi/2, A , and tan/3 = (H 2 )/(Hi), in addition to the sign of fi. The 
parameters m^, and A can be complex, and four phases appear: 6a (from Ao), 9\ (from 
the gaugino mass mi), 6*3 (from the gaugino mass m.3), and 9^ (from the fi term). 

The mSUGRA model can not explain the large negative value of sin(20!)^x s , because in 
this model, the only source of flavor violation is in the CKM matrix, which can not provide a 

8 



TABLE IV: sin(20i)0x s (left column) and B{B ± cpK^ x 10 6 (right column) at tan/3 = 10 for 



various values of the parameters |A)|, m i/2 anc l 
GeV. 



AA/1 23(32) 



. The unit for these parameters are in 



14)1 


800 


600 


400 





AA D 

^^23(32) 


m i/2 = 300 


-0.55 


9.9 


-0.57 


9.2 


-0.55 


9.1 


-0.54 


8.1 


66 - 74 


m 1 / 2 = 400 


-0.56 


9.9 


-0.53 


9.6 


-0.56 


9.2 


-0.58 


8.5 


150 - 168 


m l/2 = 500 


-0.37 


9.9 


-0.39 


9.9 


-0.42 


10.0 


-0.43 


8.1 


244 - 256 


nil /2 = 600 


-0.32 


7.6 


-0.30 


7.5 


-0.30 


7.5 


-0.05 


7.1 


270 - 304 



sufficient amount of flavor violation needed for the b — > s transition in the processes B — > <pK 
jl^ . The minimal extension of the mSUGRA has been studied to solve the large negative 
sin(20!)^x s in the context of QCDF [if]], or both large negative sin(20x)^ s and large BR 
of B -> rfK in the context of NF Q- 

We consider only non-zero (2,3) elements in A terms as a simplest extension of the 
mSUGRA model. This new piece enhances the left-right mixing of the second and third 
generation. The A terms with only non-zero (2,3) elements can be expressed as 



A U,D = A Y U > D + AA U,D , 



\U,D 



AA 



U,D 



e l<p v with 



AA 



U,D 
ij 



(19) 
= 



\U,D 



where AA ' are 3x3 complex matrices and AA^ 

unless (i, j) = (2, 3) or (3, 2). It is obvious that the mSUGRA model is recovered if AA 
0. 

It has been noticed il6j that the SUSY contribution mainly affect the Wilson coefficients 
C 8s (7 7 ) and (^85(77) an d these coefficients do not change the weak annihilation effects arising 
from the SM calculation. In our analysis, we consider all the known experimental constraints 



on the parameter space of the model, as in Ref. 



0. 



Those constraints come 



rom the 



25 



radiative B decay process B — > X s j (2.2 x 10" 4 < Br(B — > X s j) < 4.5 x 10~ 

neutron and electron electric dipole moments (d n < 6.3 x 10~ 26 e cm, d e < 0.21 x 10 _26 e cm 
22|), relic density measurements, K° - K° mixing (AM K = (3.490 ± 0.006) x 10~ 12 MeV 
271]). LEP bounds on masses of SUSY particles and the lightest Higgs {m h > 114 GeV). 

From the experimental constraints, we find that 81 ~ 22°, 83 ~ 30°, and 8^ ~ —11°. For the 

phase 8a, we set 8a = tt. 
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TABLE V: Ab~> s +~/ (left column) and A ( j > K± (right column) at tan f3 = 10. 



141 


800 


600 


400 





mi/2 = 


300 


0.011 


0.21 


0.017 


0.22 


0.02 


0.22 


0.024 


0.23 


"»i/2 = 


400 


0.017 


0.21 


0.028 


0.21 


0.027 


0.22 


0.029 


0.23 


mi/2 = 


500 


0.033 


0.18 


0.034 


0.18 


0.033 


0.19 


0.03 


0.21 


mi/2 = 


600 


0.023 


0.2 


0.023 


0.2 


0.023 


0.2 


0.02 


0.16 



TABLE VI: B(B ± -> rfK^) x 10 6 (left column) and B(B ± -► r/K ± ) x 10 6 (right column) at 
tan [3 = 10. 



141 


800 


600 


400 





mi/2 = 


300 


79.6 


3.67 


81.0 


3.76 


79.6 


3.66 


79.0 


3.62 


mi/2 = 


400 


78.2 


4.4 


83.0 


3.85 


79.0 


3.69 


81.0 


3.69 


^1/2 = 


500 


84.8 


3.90 


83.7 


3.85 


81.0 


3.71 


77.0 


3.50 


mi/2 = 


600 


73.0 


3.26 


71.0 


3.18 


70.0 


3.10 


70.0 


3.00 



Now we consider the case with non-zero AAg and AAg for tan/3 = 10 in our calculation 
(Similar results can be obtained for other values of tan/5). All the other elements in AA U,D 



are set to be zero. We also implement 



AAg 



AAg and (pg 3 ^ <pg. In this case, in order 



to find solutions, we do not need to use the additional strong phase: i.e., 5' = 0. The value of 
TO1/2 varies from 300 GeV to 600 GeV, and the value of \Aq\ varies from to 800 GeV. Even 
though the value of m is not explicitly shown, it is chosen for different and A such 
that the relic density constraint is satisfied. (We satisfy the relic density constraint using the 
stau-neutralino co-annihilation channel The value of m increases as mu2 increases. 



The value of 



/\A D 



increases as m\/2 does. The phases 0^ and (j)^ 2 are approximately 
-40° to -15° and 165° to 180°, respectively. 

Table IV shows the BRs for B — > <f>K (right column) and sin(20i)0x s (left column), 
calculated for various values of the parameters mi/2 and \A \. The value of sin(20x)^ s 
can be large and negative. If we use the additional strong phase 5' ^ 0, the magnitude 
of sm{2(j)i) ( i >Ks can be even larger. For example, if 5' = —30°, then sin(20x)^ s = —0.65 
for mi/2 = 400 GeV and Aq = —800 GeV. In Table V, we show the CP rate asymmetries 
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TABLE VII: A^x+ (left column) and A v k+ (right column) at tan (3 = 10. 



\M 


800 


600 


400 





m 1/2 = 300 


0.0 


-0.24 


0.0 


-0.23 


0.0 


-0.23 


-0.001 


-0.22 


m \/2 = 400 


-0.026 


-0.31 


-0.003 


-0.22 


-0.004 


-0.24 


-0.003 


-0.23 


m \/2 = 500 


-0.005 


-0.21 


-0.005 


-0.21 


-0.008 


-0.22 


-0.002 


-0.22 


m 1/2 = 600 


-0.0004 


-0.22 


-0.0006 


-0.22 


-0.001 


-0.22 


-0.009 


-0.195 



Ab^ s +-y (left) and A^k± (right). Since A 23 contributes to b — > 57, the CP asymmetry gets 
generated, but the asymmetry is small. So far we have assumed that AA^ 3 32 = 0. But if 
we use AA 2 / 3 32 7^ and AA 2 : 3 3 2 = 0, the value of sin(20)^x s is mostly positive. 

The BRs for — > i]'K ± and — > r]K ± are shown in Table VI and the CP asymmetries 
are shown in Table VII. They are consistent with the experimental data. The value of 
sin(20i)^ is 0.6-0.7 for all these scenarios. The BRs are: B(B° -> rj'K ) ~ (63-73) x 10" 6 
and B(B° -> rjK°) ~ (2.25 - 3.37) x 10~ 6 . The CP asymmetries are: A v > K o ~ A^ K + and 
A v k° ~ A r) K+- We see that the CP asymmetries are also in good agreement with the 
experimental values shown in Table I. As a final comment, we note that in the case of 
the large Xa,h effects with Xmm = 7.5, it is possible to raise the BR for B + — > <f)K + to 
(8 - 9) x 10~ 6 . However, in that case, A C p(B + -> <f)K + ) ~ -0.2 and sin(2^i)^, > -0.34. 

In conclusion, we have analyzed the CP asymmetries of B — > <pK and B — > modes 
in the QCDF framework. The phenomenological parameters X H and Xa arising from end- 
point divergences in the hard spectator scattering and weak annihilation contributions are 
determined by the global analysis for twelve B — > PP and VP decay modes, such as 
B — > 7T7T, ttK, pn, pK, etc, but excluding the modes whose (dominant) internal quark- 
level process is 6 — > sss. We found that it is possible to explain large negative sm(2(fri) ( f )K3 
simultaneously with the large BRs for B ±< - ' ) — > rj'K^ ^ in the context of supersymmetry. 
We use R-parity conserving SUGRA models and constrain the parameter space with the 
observed dark matter relic density along with other experimental constraints. We also show 
that the R-parity violating SUSY models can provide solutions. Our results are consistent 
with other available data on CP asymmetry for different B — > <j)K, B — > rj'K, and B — > r]K 
modes. 
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